Size of direct CP violation in singly Cabibbo-suppressed D decays 
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The first experimental evidence for direct CP violation in charm-quark decays has recently been 
presented by the LHCb collaboration in the difference between the D — »> K + K~ and D — >• 7r + 7r - 
time-integrated CP asymmetries. We estimate the size of the effects that can be expected within 
the Standard Model and find that at leading order in l/m c they are an order of magnitude smaller. 
However, tree-level annihilation type amplitudes are known to be large experimentally. This implies 
that certain formally l/ra c -suppressed penguin amplitudes could plausibly account for the LHCb 
measurement. Simultaneously, the flavor-breaking parts of these amplitudes could explain the large 
difference between the D — > K + K~ and D — > tt + tt~ decay rates. 



The LHCb collaboration has recently announced 
a measurement of the difference between the time- 
integrated CP asymmetries in two singly Cabibbo- 
supressed (SCS) D-meson decay modes, Acp(D — >• 
K+K-)-A C p(D 7T+7T-) = (-0.82±0.21±0.11)% ffl. 
In the standard model (SM), CP violation (CPV) in D 
decays is commonly expected to be very small. It has 
often been stated that "a measurement of CP violation 
in D decays would be a signal of new physics". If this 
statement holds under further scrutiny, the LHCb mea- 
surement may well be the first evidence for new physics 
(NP) at the LHC. The measurement can be identified, to 
excellent approximation, with the difference between the 
two direct CP asymmetries 

AAcP = Af +K --A% v - (1) 

(the contribution of indirect CPV, due to a small differ- 
ence in the cut on the D° proper decay time in the two 
decays, is constrained to be < 0.03%). The new world 
average for the direct CP asymmetry difference becomes 
AA C p (-0.67 ± 0.16)% 0, after the inclusion of the 
most recent CDF, BABAR and Belle measurements [3]- 
[5]. We are therefore motivated to examine the size of 
direct CPV in the SM. 

The role of the third generation in D meson decays 
and D° — D° mixing is suppressed by small Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements. Thus, 
all observables are CP conserving to first approxima- 
tion. More precisely, in the SM CPV in D° — D° mix- 
ing arises at 0(V ch V ub /V cs V us ) ~ 1(T 3 . Direct CPV 
in SCS decays is further parametrically suppressed to 
0([V c bVub /V cs V us \a s i / 7r) ~ 10 -4 , since it comes from the 
interference of the tree and penguin amplitudes. Does 
this naive scaling estimate for AAcp mean that NP has 
been discovered? It is interesting that despite severe con- 
straints from other FCNC processes, D—D mixing in par- 
ticular, it is possible to attribute AAcp entirely to NP 
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[6]. As far as the SM predictions are concerned, in [7] it 
was argued that large CP asymmetries can be expected 
in SCS D decays, while in [8 small CP asymmetries were 
obtained. Since these early attempts, much more exper- 
imental information on D decays has become available. 
In this paper we address the following question: are there 
any indications in the data that AAcp is due to SM con- 
tributions that are enhanced well above the naive scaling 
estimate? 

We will show that there is a possible dynamical ex- 
planation for AAcp m the SM. A requirement for it to 
be viable is that the two matrix elements of the "tree" 
operators in the effective weak Hamiltonian, 

(K+K-\(dd) V - A (uc) V - A \D°), (2) 

and 

(K + K-\(ss) V - A (uc) V - A \D°), (3) 

are of similar magnitude. (For simplicity we show only 
one of the two possible color structures in Q and 
The analogous requirement applies to the matrix ele- 
ments obtained via the replacements K — >• ir and s «-»• d. 
The matrix element ([2| only contains contributions in 
which the d and d fermion fields are contracted in the lo- 
cal operator. The matrix element ([3| receives contracted 
as well as uncontracted contributions, i.e. the contribu- 
tions where the s and s fermion fields are contracted with 
the final state kaons. Our requirement is therefore equiv- 
alent to saying that the contracted contributions need to 
be enhanced over the uncontracted ones. 

We will give two indications that this is the case. First, 
assuming that the two matrix elements ([2| and ([3| dif- 
fer by 0(JkI fir — 1) ~ 0.2, a quantity identified with 
nominal SU(3) breaking, it is possible to simultaneously 
explain both the large ratio Br(D — » K Jr K~)/Br(D — » 
7T + 7r - ) 3.3 and the enhancement of AAcp- Second, 
using a one-gluon exchange approximation, we identify 
formally power-suppressed contracted contributions in 
Q and Q that indeed show dynamical enhacement. It 
is important to note that the enhancement of contracted 
over uncontracted contributions in Q and Q is in prin- 
ciple measurable on the lattice [9]. Therefore, the SM 
explanation of large AAcp can be definitively checked 
in the future. 
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We begin with a discussion of enhanced contractions 
in the one-gluon exchange approximation. The essential 
ingredients are: (i) 1/N C counting; (ii) D branching ra- 
tio data which shows that certain formally 1 / m c power- 
suppressed amplitudes are of same order as their leading 
(l/m c )° counterparts; (hi) translation of this breakdown 
of the l/m c expansion to the penguin contraction am- 
plitudes, in the approximation of a hard gluon exchange; 
(iv) a perturbative estimate of the related "effective Wil- 
son coefficients". 

The starting point in estimating A^\ K - and is 
the effective weak AC = 1 Hamiltonian i^ e ff, evaluated 
at scales \i ~ m c . It is obtained by integrating out the W 
boson and running at NLO down to \i ~ m^, where the 
6-quark is integrated out using NLO matching. Finally, 
running down to /i ~ m c yields [TO] 

#e A ff c=1 = % [ E k v «p ( c iQi + cm 

p=d,s 



(4) 



H.c. 



i=3 



The "tree" operators are Q\ = (pc)v-A(up)v-A, Q\ — 
(PaCp)v-A(upp a )v-A with summation over color indices 
a,/3 understood. The penguin operators are Qs y ^ = 



and Qsg = 



8tt' 



mcUcr^vil + 75) G^c, where the sum- 



mation is over the light quark species, q = u,d,s. The 
Wilson coefficients Ci(ja) are evaluated at \i ~ m c . We 
find d = 1.21, C 2 = -0.41, C 3 = 0.02, C 4 = -0.04, 
C 5 = 0.01, C 6 = -0.05, C 8p = -0.06. 

The CP-conjugate decay amplitudes for SCS D — >> / 
and D —> f decays in the SM can be written as 



Af(D -> f) = Aj'[l + r /e * ( ^- 7) ] 



(5) 



where A J is the dominant tree- level amplitude, taken to 
be real and positive by convention, and rj = A^/A^ 
is the relative magnitude of the subleading "penguin" 
amplitude, Aj , which carries the CKM weak phase 

7 = (67.313 5 )° i 11 .- The penguin amplitude is CKM 
suppressed by 0(V c bV u b/V cs V us ) compared to the tree- 
level amplitude and carries the relative strong phase Sf. 
The direct CP asymmetry is then 



\Aj? 



l%l 2 



= 2rf sin 7 sin (5/, 



(6) 



where we only kept the term linear in 7*/. 

From the measured D — >• / branching ratios we can 
obtain information on the sizes of the tree-level ampli- 
tudes A J. It is important to note that the data on SCS 
or Cabibbo favored (CF) decays implies that formally 
l/m c -suppressed tree- level amplitudes are of the same 
order as the leading power contributions. For instance, 



the tree level CF decay D° 
diagrammatic notation, is 



K 7r + , written in SU(3)f 



A T 



K+7T- 



V* s V u d(T K7 r + E K 7t), 



(7) 



and similarly for the tree-level amplitudes in D° — >• 
K + K~, 7T + 7r - with appropriately modified CKM ele- 
ments. Here Tf are the "color- allowed VF-emission" am- 
plitudes, and Ef are the "W-exchange" annihilation am- 
plitudes p~2j [13] . Equating the magnitudes of the total 
tree amplitudes to the measured ones, one has Aj c _ 7r+ ~ 
2.5 keV ~ A^./X ~ A T K+K _/\, with A = \V US \ =0.22. 

At leading power and in naive factorization Tf take 
the familiar form, e.g., T n7r ex Fjj^^f^^ where Fjj^^ 
is the form factor and f n the decay constant. The an- 
nihilation amplitudes Ef are formally l/m c suppressed 
power corrections. However, SU(3) fits to the data yield 
(/ = KK, tttt, Kit) HUTS] ' 



Ef ~ T f and E KK ~ E n7T ~ E K7r . 



(8) 



The first relation signals the breakdown of the l/m c ex- 
pansion. 

The QCD penguin amplitudes are 

A P k+ k- = -V: b V ub P KK , < +7r _ = -V c * b V ub P nn , (9) 

yielding weak phases —7 and 7r— 7 relative to the D — >• tttt 
and D — > KK tree amplitudes, respectively, with sin 7 ~ 
0.9. This implies that the CP asymmetries in D — >• tttt 
and D — >• have opposite signs. 

First let us check whether the naive CKM scaling esti- 
mate rf ~ O(10~ 4 ) is reproduced in the formal m c — > 00 
limit. Rough estimates for these leading-power penguin 
amplitudes are obtained using QCD factorization at NLO 
in a s [6j[16l[T7], i.e. in naive factorization plus the 0(a s ) 
corrections coming from down- and strange-quark loop 
penguin contractions, vertex corrections, and hard spec- 
tator interactions. For the penguin to measured tree am- 
plitude ratios 



r f 



A^ (leading power) 



Aj(exp) 



(10) 



we obtain r^ +K _ w (0.01 - 0.02)%, r; +7r ~ 
(0.015 — 0.028)%, for a renormalization scale in the range 
fi G [lGeV, m^], and down and strange quark masses 
in the penguin contraction loops varied in the ranges 
m^,m s ~ 0.1 — 0.4 GeV. This is consistent with, yet 
slightly bigger than, the naive expectation that rf ~ 
0([V cb V ub /V cs V us ]a s /7r) ~ 10- 4 , based only on CKM 
factors and a s scaling. Assuming 0(1) strong phases 5/, 
the direct CP asymmetries satisfy A^ 1T ~ 2r/. The U- 
spin relation A^ = —Aj^ r K then leads to AAcp ~ 4r/, 
or AA C p = 0(0.05% - 0.1%), an order of magnitude 
smaller than the measurement. However, already at lead- 
ing power our estimates for the penguin amplitudes are 
dominated by the penguin contractions. 
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FIG. 1: Annihilation topologies. 




FIG. 2: Penguin contraction topologies. 



Next, we turn to the QCD penguin power corrections. 
For concreteness we focus on the annihilation topology 
amplitudes, which in the tree amplitudes are seen to 
be of the same size as the leading power contributions, 
Tf ~ Ef. We consider two examples of l/m c penguin 
amplitudes for the final states / = K + K~ and 7r + 7r - , 

P fil =^ G 6 </|- 2 (uu) s+ P ® A (uc) S -p\D), (11) 



and 



P />2 =-j= (C 4 + C 6 )(/| (q a qp) V ±A ® A (u c a ) v _ A \D), 

(12) 

corresponding to insertions of the QCD penguin opera- 
tors in the left and right annihilation topology diagrams 
of Fig. [I] respectively. Summation over q = u, d, s is un- 
derstood. The annihilation product J\ A J2 means that 
J2 destroys the D meson, and J\ creates the final state. 
We will estimate the non-perturbative matrix elements 
in (11) and (12) from the sizes of Ef in ([7]), see below. 

The penguin amplitudes Pf^ and Pf^ also receive con- 
tributions from penguin contractions of the tree operator 
Qi, shown in Fig. [2] (left) and Fig. [2] (right), respectively. 
These contributions cancel the log \i scale dependence of 
the Wilson coefficients in ( jTT| ) and ( [l2| ). In the partonic 
picture the cancelation in Pf^ (Pf^^jis associated with 
the exchange of a single gluon between the s and d-quark 
loops and the spectator quark (outgoing qq pair), but 
an arbitrary number of partons otherwise exchanged be- 
tween the external legs. This leads to effective coefficients 
for 64,6 that depend on the gluon's virtuality q 2 , 



Cq (4) (aO 



rl 
6 



log 



1 



2tt 



8 



-G 



( m l m d q 2 \ 



(13) 



where G(si,S2,x) = (G(si,x) + G(s2,x))/2 is an in- 
frared finite quantity. At one loop, G(s,x) is defined, 
e.g., in [17]. 

In order to obtain rough values for the penguin con- 
traction amplitudes we make two approximations: (i) we 
use the partonic quantity G as an estimator for the un- 
derlying hadronic interactions, for instance, final state 
interactions, (ii) we take G at a fixed value of q 2 = q 2 . 
With these approximations we can replace C^q with C%% 
in (11) and (|12|). Order of magnitude estimates of the 



annihilation- type matrix elements in (11), (12) are ob- 
tained by appealing to the 'W-exchange" annihilation 
amplitudes Ef in ([7]). For instance, for D° — >• i^ + 7r~, we 
have 



<g) A (upc a ) 



V-A 



\D) 



(14) 



neglecting a smaller contribution proportional to C2 ~ 
-Ci/3. Note that the E f enter at 0(1/N C ). 

We are now able to study the ratios 77 in Eq. (10) 
beyond leading power, 



.PC 



A j (power correction) 



Aj(exp) 



, 7T + 7T 



(15) 

by assuming N c counting, and ignoring differences due to 
the different chirality structures, which should be small 
in D — )> PP decays. We thus equate the operator matrix 
element fL4|) with the matrix elements of the operators in 
(11) and (|12|), taking into account that Pf^ is 0(1) and 
P/ ? 2 is 0(T/N c ). We use the experimentally supported 
relations E f ~ T f ~ Aj(exp), cf. Eq. Finally, to 

arrive at numerical estimates we set Ef = Aj(exp) and 
use the experimental observation 



Ajf^Oc to obtain 



AT f 

A f,2 



2N c \V cb V ub Cf\/(C 1 sin# c ), 



(16) 



2\V cb V ub (Cf + C 6 eff )|/(C lS in^ c ). (17) 



The above estimates for 77 depend on q 2 through the 
G function. Taking q = m c , varying \i G [lGeV, m^] 
and setting m s = 0.3 GeV, rrid = 0.1 GeV we obtain 
77,1 G [0.04,0.08]% and 77,2 G [0.03,0.04]%. We see 
that the individual power corrections which we have con- 
sidered could be enhanced by a factor of a few relative 
to leading power. For example, for a renormalization 
scale jjL = 1 GeV and taking AAcp ~ 4r/, we find 
AA C p ~ 0.3%(0.2%) for P fy i(P fy2 ). Assuming smaller 
constituent quark masses could give even larger effects. 

The enhancement of the penguin amplitude can be 
traced to (i) the breakdown of the l/m c expansion, via 
the comparison of the matrix elements in P/,2 an d 
Ekk, and (ii) the penguin contraction contributions, 
which dominate in (13) via the G function. 
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There are other contributions to the penguin ampli- 
tudes at subleading order in l/m c which may be at least 
as large as those we have discussed. An example is pen- 
guin contractions which correspond to emission of two or 
more gluons from the s and d-quark loops in the partonic 
picture (these were discussed for B — >> Kir in [18]). 

We now address the other indication for enhanced pen- 
guin contractions. Assume that this enhancement is of 
0(1/ e), where e - 0(20%) is the typical size of SU(3) 
breaking. Then 

A(D° -> 7T+7T-) =V cd V: d (T + E + P breaU ) 
- V* b V ub P, 

A(D° -> K+K-) =V CS V; S (T + E- P break ) 

-v; b v ub p, 

where we only write down the leading contributions in e 
counting for each product of CKM elements. The O(e ) 
amplitudes T and E equal Tk-k an d Ektt in in the 
U-spin limit. The U-spin symmetric penguin amplitude 

p = W\ E CiiQf + QDlD ), (20) 

2=1,2 

where / = K + K~ or 7r + 7r - , contains contracted contri- 
butions and is thus 0(l/e). The U-spin breaking differ- 
ence of the penguin contractions, 

Pbreak = \ £</| £ C ^ ~ Q\W) , (21) 

is part of the 0(e - 1/e) ~ O(e ) amplitude. In partic- 
ular, it gives 0(1) corrections to the A (D° — > 7r + 7r~) 
and A(D° — >• K + K~) amplitudes, explaining the large 
difference in the branching ratios (see also [19]). 

At the same time P ~ Pbreak/e is of the correct size to 
give 



VcbVub 




P 






T + E ± Pbreak 


\v cb v ub \ 


1 


-0.15%, 


\v cs v us \ 


2e 



(22) 



for e — 0.2. For 0(1) strong phases this gives an esti- 
mate |A^4cp| ~ 4r/ ~ 0.6% that is strikingly close to 
the world average [2] (for further details see [20]). This 
shows that for nominal SU(3) breaking of (9(0.2), the 
enhancement of the contracted amplitudes can simulta- 
neously explain the large difference of Br(D K + K~) 
and Br(D — >> 7r + 7r _ ), and the SM enhancement of AAcp 
over the naive scaling estimate. 

Given the measured value of AAcp, can we predict the 
orders of magnitude of direct CPV in other SCS D de- 
cays? For processes that only differ from D° — >> K + K~ 
and D° — )> 7r + 7r _ by a change in the flavor of the spec- 
tator, we would generically expect direct CP asymme- 
tries of same order, both in the SM and in NP exten- 
sions that modify the QCD penguin operators. Exam- 
ples are D + K+K° and £>+ -> tt+K , which contain 



SM contributions from the penguin contractions of the 
type shown in Fig. [2] (left) or but not Fig. [2] (right) 
or .P/,2- Interestingly, for D — >• KsKs the rate is due 
to the difference of s and d quark penguin contractions 
of the latter type, as well as the difference of exchange 
graphs E with outgoing s and d quark pairs, while the 
SM penguin amplitude is due to the two sums. Thus, a 
large direct CP asymmetry is possible in the SM for 0(1) 
strong phases, Af sKs - \V cb V ub /V cs V us \2/e - 0.6%. A 
study of QCD penguin amplitudes in SCS decays to 
vector and pseudoscalar pairs is also of interest. 

It is of course not excluded that AAcp is (at least 
partly) due to NP. In fact, it is possible to prove the 
presence of NP in AI = 3/2 decays. For example, the 
D + — ^ 7T + 7r° amplitude does not receive a QCD penguin 
contribution. It must vanish because 7r + 7r° is an I = 2 
final state which cannot be reached from the initial state 
via the AI = 1/2 QCD penguin operators. It does, how- 
ever, receive an electroweak penguin contribution which 
is a/ as suppressed in the SM. Thus, enhanced direct 
CPV in D + — )> 7T + 7r°, e.g., at the level of the observed 
AA CP , would be a signal for isospin violating NP [2T] , 

In conclusion, we have shown that it is plausible that 
the SM accounts for the measured value of AAcp- At 
leading order in l/m c we obtain AAcp = 0(0.05% — 
0.1%), an order of magnitude smaller than observed. 
However, we find that significant enhancements are pos- 
sible for certain penguin contraction power corrections 
whose matrix elements can be estimated using informa- 
tion obtained from global SU (3) fits to the measured D 
decay branching ratios, utilizing N c counting for guid- 
ance. We reach the conclusion that the ratio of QCD 
penguin to tree amplitudes can be rf ~ 0.1%. Us- 
ing the /7-spin relation A dl + _ — — Aj^ r +K - one has 
AAcp — 2A dir — 4=rf sin S sin 7. For a large strong phase 
S - 0(1) this yields AA C p ~ 0.4%. This result is sub- 
ject to large uncertainties due to (i) the extraction of the 
tree- level annihilation amplitudes Ef from data, (ii) ne- 
glected contributions in Ef, (hi) the use of N c counting 
for penguin operator matrix elements, (iv) the modeling 
of Qi penguin contraction matrix elements, and (v) ad- 
ditional penguin contractions, not associated with log \i 
cancellations, that have not been included in our esti- 
mates. A cumulative uncertainty of a factor of a few is 
reasonable. 

Further support for the above SM explanation, where 
AAcp is due to the enhancement of contracted Qi op- 
erator contributions, is provided by the large difference 
Br(D K+K~) ~ 2.8Br(D tt+tt-). It can be 
attributed to the SU(3) breaking part of the penguin 
contractions, if these are enhanced. This fixes the size 
of the penguin contractions and leads to the prediction 
Tf ~ 0.2%, in agreement with our semi-perturbative es- 
timates, and most importantly, in agreement with the 
measured value of AAcp for strong phases of order (9(1). 
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